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Abstract 

We compute by Density Functional Theory-Non Equilibrium Green Functions Formalism (DFT- 
NEGFF) the conductance of bulk Si along different crystallographic directions. We find a ceiling 
value for the intrinsic mobility of bulk silicon of 8.4 • 10 6 cm 2 /Vs. We suggest that this result is 
related to the lowest effective mass of the <001> direction. 

Silicon is the most widely used material in current electronic devices. Due to downscaling, much work 
has been done on nanostructures such as nanowires. Their structural and electrical properties have 
been studied experimentally and theoretically pQ and devices have been made with them [2]. Since 
these new structures are to be used in electronic devices, an understanding of their electronic transport 
behaviour is essential. 

Carrier mobility is an important transport property of a semiconductor. It reflects the capacity of the 
charge to move in the material and to respond to field changes [3]. 

Due to the challenges of current ah initio computations, transport properties have been studied by this 
method only in very low dimensional systems [H [6] . The most common way to examine transport 
and mobility in bulk silicon has been to solve the Boltzmann equation using MonteCarlo methodology 
with electron-phonon interaction as a scattering mechanism [HJ, O [TU1 lllj . To our knowledge, no 
theoretical approach based on fully atomistic ah initio methods has been used to study bulk behavior. 
Materials in devices are normally used in an anisotropic way: they are grown along a certain axis 
and one direction of transport becomes predominant, for example, silicon nanowires (SiNWs). The 
transport properties are expected to be sensitive to the crystal orientation, due to the anisotropy of 
band structure of bulk silicon. Thus, a detailed understanding of the electrical properties as a result 
of the different growth conditions is required in order to assess the impact of the crystal structure 
on transport properties and to determine which orientation is best suited to the engineering of high 
performance thick SiNW-based devices. 

In this Letter, we present a theoretical study of the electron transport properties of bulk silicon as an 
approximation to the transport in thick nanodevices, where the bulk properties are still predominant 
over confinement effects [T2J. For the sake of completeness, the most significant growth directions 
of SiNWs, namely cubic diamond phase oriented along <001>, <110> and <111> axis, have been 
considered. We have evaluated the intrinsic electron mobility due to elastic scattering of the electrons 
inside the crystal. We compare the values with current experimental results from the literature and 
with the corresponding results for small SiNWs, finding a ceiling of the value of bulk silicon mobil- 
ity depending on the crystallographic direction which is lower than the reported value for SiNWs 
[131 HH, [15] . The authors suggest that surface states may be responsible for the measured conduction. 
Here the electronic transport is studied inside an open-boundary bulk system described by first- 
principles field within Non Equilibrium Green Functions Formalism (NEGFF) |16|. H7] . We performed 
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Density Functional Theory (DFT) calculations with a numerical atomic orbital basis set [18]. We 
used a double-C polarized basis set, with the generalized gradient approximation [19] for the exchange- 
correlation functional in the PBEsol parametrization [20]. Real space mesh cut-off and k-grid values 
were set to converge total energy values within 0.1 meV and maximum absolute forces over atoms are 
smaller than 0.004 eV/Ang in all cases. 

The electronic transport properties were studied with Non Equilibrium Green Functions Formalism, 
within the Keldysh formalism |21j . based on DFT as implemented in TranSIESTA code [22 . It com- 
putes the electronic transport properties due to elastic scatter in a ballistic regime of a nanoscale 
system confined in a central scattering region (C) in contact with two semi-infinite left (L) and right 
(R) electrode regions at different electrochemical potentials. The bias voltage Vb between them corre- 
sponds to the electrochemical potential difference between the left (up) and the right (//r) electrodes, 
eVb = up — ur. In our case, both electrodes and the contact region were made of Si in the same 
crystallographic orientation. 

The density matrix from the incoming scattering states of the left (1) and right (r) electrode is 
D (x,y) =^2ipi(x)^*(y)n F (ei - Hl) +J2A(x)^*(y) n F(e r - m), 

l r 

where the scattering states starting in the left electrode are 

^{x) = + / dyG(xy)V L (y)^(y). 

These are generated from the unperturbed incoming states of the uncoupled, semi-infinite electrode 
ip®, using the retarded Green's function G of the coupled system [22J and the interaction between the 
left electrode and the contact region Vp. 

The current I through the contact region was calculated using the Landauer-Buttiker formula [23], 

/oo 
de[n F {e - u L ) - n F (e - u R )]T(e, V b ), 
-oo 

where Go = 2e 2 /h is the unit of quantum conductance, np is the electrode distribution function and 
T(e, Vb) is the transmission probability of electrons incident at an energy e through the device under 
the potential bias Vb- 

We have computed the I(V) characteristic (Fig. [TJ in each direction and evaluated the intrinsic con- 
ductance. 

We propose to calculate the corresponding electron mobility fx of the electrons along the different di- 
rections of bulk Si. Assuming an intrinsic concentration of carriers n corresponding to a temperature 
of 300K, the corresponding mobility \i = a /en (where a is the associated conductivity computed as 
the conductance across the corresponding unit cell in each direction) is shown in Fig. [2] 

The mobility in bulk semiconductors as we go down from room temperature increases at first due to 
the suppression of phonon scattering [16] . But it does not increase any further once the phonon scat- 
tering is small enough so that impurity scattering becomes the dominant mechanism. With undoped 
samples we can obtain higher mobilities, because we have suppressed inelastic scattering centers (Fig. 

In agreement with these considerations our computed mobility is greater than experimental results 
and thus constitutes a ceiling value for them. TranSIESTA code only takes into account the transport 
due to elastic scatter inside the crystal, without considering the contribution of inelastic scatter due to 
phonons and impurities. Despite current advances [25], up to now there is no reliable way to compute 
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Figure 1: Comparison of the I(V) characteristic (I in this case is the current per unit of transversal 
area) of the three directions of cubic diamond <001>, <110> and <111>. 
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Figure 2: Extract from ref. [3]. Experimental mobility of electrons in Si as a function of temperature 
with different impurity concentrations (point symbols). The dash-dot line indicates the pure lattice 
mobility obtained with Monte-Carlo simulation [9]. Our computed mobilities (continuous lines: red 
for <001>, blue for <110> and green for <111>) are compared with experimental values for Si bulk 
(point symbols [3] and square symbol |24j ) and experimental [14J (diamond symbol) and theoretical 
(with NEGFF) [15] (dot-dot line) values of thin SiNWs. 
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inelastic scattering in this framework. On the other hand, our calculated conductivities are lower than 
experimental ones in SiNWs because in bulk the surface states are not taken into account, and they 
seem to play an important role in transport properties in very thin SiNW [Tl 1141 [T5] , Thus, our result 
becomes a threshold to distinguish bulk conductivity from conductivity affected by nanoscale effects. 
Electrons in silicon occupy several valleys which have different orientations. The effective mass associ- 
ated with these valleys gives the response of the electrons to an external field; the electrons with lower 
effective mass are faster, and the direction of these valleys is the preferred direction of the electron 
transport. The results are similar to the experiments |26j to increase the mobility due to an applica- 
tion of a uniaxial stress along a direction such that the valleys presenting the lower effective mass in 
the direction of the field have less energy than the other valleys. Electrons in equilibrium will have 
greater tendency to occupy these "faster valleys", and the overall mobility will increase. 
We have calculated the effective mass (see Table [T]), with respect to the free electron mass mo, around 
the conduction band minimum (CBM) in each direction using the expression AE c (k) = h 2 k 2 /2m*, 
where m* is the electron effective mass around the CBM. We fit a quadratic of the form y = mx 2 , then 
we solve m = fr 2 /2m*. For bulk Si in <001> direction, we calculate m* to be 0.884 mo- This value 
corresponds to the longitudinal effective mass of electrons in bulk silicon and it is close to the reported 
experimental values [27 1 . In order to check the calculations we computed the electron transversal ef- 
fective mass and the light hole effective mass, and the values were 0.194 mo and 0.163 mo, respectively. 
They are close to the reported experimental values [28]. 





<001> 


<110> 


<111> 


electron effective mass 
mobility (cm 2 /Vs) 


0.884 m 
8.4 • 10 6 


1.198 m 
6.04 • 10 6 


1.575 mo 
3.99 • 10 6 



Table 1: Electron effective mass and mobility for bulk Si oriented in <001>, <110> and <111> 
direction. 

In conclusion, we have performed a theoretical study of the mobility of bulk Si based on the description 
with NEGFF of the electronic transport. We have evaluated different crystallographic directions and 
found the <001> the best for electrons in bulk Si due to the absolute minimum of the lowest conduc- 
tion bands in the band diagram and the lower corresponding effective mass. We have performed the 
calculations inside the framework of Landauer-Buttiker formalism and found an intrinsic conductance 
due to elastic scattering of electrons inside silicon bulk, and we can associate a mobility value to it. 
In this sense, our computed values are a ceiling (i.e. the maximum possible mobility that the material 
could achieve). With the contribution of inelastic scatter due to phonons and impurities, the mobility 
is reduced. Moreover, our result becomes a threshold between bulk conductivity and conductivity 
affected by nanoscale effects. While the Non Equilibrium Green Functions Formalism is currently 
used to study transport in very small nanostructures, we have used it in huge configurations (bulk), 
showing that it is a very powerful technique for evaluating properties of many-body systems. 
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